Tracer diffusion in a dislocated lamellar system 
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Many lamellar systems exhibit strongly anisotropic diffusion. When the diffusion across the 
lamellae is slow, an alternative mechanism for transverse transport becomes important. A tracer 
particle can propagate in the direction normal to the lamellae, never leaving a particular layer, by 
going around a screw dislocation. Given the density of positive and negative screw dislocations, 
we calculate the statistical properties of the transverse transport. When either positive or negative 
dislocations are in excess, the tracer moves ballistically normally to the layers with the mean square 
of the displacement growing like the square of time T 2 . When the average dislocation charge is 
zero, the mean square of the normal displacement grows like T log T for large times. To obtain this 
result, the trajectory of the tracer must be smoothed over distances of order of the dislocation core 
size. 
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Diffusion in layered systems is often strongly 
anisotropic. The mechanisms and the manifestations 
of the anisotropy vary. In many such systems diffu- 
sion across is significantly slower than along the lamel- 
lae. For example, enhanced creep resistance of lamellar 
alloys, such as industrial TiAl, is probably due to the 
high barrier for the dislocations crossing from one layer 
to the next iQ §]. In lamellar phases of diblock copoly- 
mers ||, tracer diffusion along the lamellae can be up 
to forty times faster than across [Q. The fact that wa- 
ter diffusion in lamellar phases of phospholipid bilayers 
Q is strongly anisotropic may be relevant to attempts 
to use multilamellar vesicles for drug delivery |^|. An- 
other example of anisotropic diffusion is the kinetics of 
electroactive probes in lyotropic liquid crystals 0. 

When lateral diffusion is much faster than transverse 
diffusion, the tracer can still be transported quickly in the 
direction normal to the lamellae if screw dislocations are 
present in the system. A screw dislocation is constructed 
by cutting a perfect layered structure with a half-plane 
normal to the layers, shifting the two sides of the cut 
with respect to each other in the direction normal to the 
layers by a distance equal to the layer spacing, and finally 
gluing the cut. Screw dislocations are indeed common 
in a variety of layered systems ||, g, [To). A summary 
of various dislocation properties in lamellar systems is 
presented in Ref. |p"lf . 

When a tracer particle confined to a particular layer 
encircles a screw dislocation, it finds itself in one layer 
higher (or lower). A tracer particle can then reach any 
point in the system while remaining confined to a layer. 
The trajectory of the tracer projected onto a plane paral- 
lel to the layers is a two-dimensional random walk. Upon 
completing a closed 2D trajectory, our random walker 
moves up or down the number of layers equal to the dis- 
location charge enclosed by the trajectory. We obtain an 
expression for this quantity by noting that when a sin- 



gle screw dislocation is present in the system, the layer 
number of the walker is the winding angle around the 
dislocation divided by 2ir. Consider now the sum of the 
winding angles around all the dislocations (the signs of 
the individual winding numbers are determined by the 
charge of the dislocations). This quantity changes con- 
tinuously. The change in this quantity along an open 
trajectory depends on the shape of the sample due to 
the contributions of the winding numbers around distant 
dislocations. However, when the walker returns to the 
origin, the change in the total winding number is the 
dislocation charge enclosed by the trajectory. Thus we 
identify the total winding number divided by 2ir with the 
layer number n(t) which for closed trajectories coincides 
with the normal displacement of the walker. 

In this letter we study the diffusion of a tracer particle 
confined to the lamellae. The tracer starts at the origin of 
layer n = at time t = and explores the x-y plane with 
diffusivity D. Let there also be a random distribution of 
positive and negative screw dislocations with densities /+ 
and /_ respectively. Our goal is to determine the nature 
of the transport normal to the layers by predicting the 
result of the following experiment. If some amount of the 
tracer material is placed at the origin at time t = 0, what 
is the density of the resulting cloud of tracer particles as 
a function of time? 

To accomplish this task we look at paths which start at 
the origin O (see Fig. [l]) at time t = and arrive to point 
E located a distance R from the origin at time T. We 
seek to define the layer number n(R, T). Since the layer 
number change is only well defined for closed trajectories 
we fix n(R, T) by completing the path r(t) with a straight 
segment OE connecting this point to the origin. We can 
then define n(R,T) to be the total dislocation charge 
enclosed by this trajectory. 

We now seek to average powers of n over positions of 
dislocations and random walks which end at r = R at 
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time t = T . We denote the average over positions of dis- 
locations with an overbar and average over random walks 
by angular brackets ( • } . Because changing the shape of 
the completing segment adds a constant to n(R,T), its 
average {n(R, T) ) has no physical meaning. However, its 

standard deviation a(R,T) — ^ n 2 ^ — {n) 2 is indepen- 
dent of the shape of the completing segment. It gives the 
size of the spreading tracer cloud at time T and distance 
R from the origin. 

We identified two qualitatively different cases. 

When /+ ^ /-, we are able to obtain a(R, T) thereby 
predicting the tracer density profile within the spreading 
cloud. This is possible because the average total dislo- 
cation charge within a closed trajectory which is propor- 
tional to its signed area. In this case we find that the 
vertical size of the tracer cloud grows linearly in time, 
i.e. there is superdiffusion across the layers. Moreover, 
the spreading cloud acquires a biconcave shape, since 
a(R,T) oc D 2 T 2 + 2R 2 DT (here brackets denote aver- 
aging over random walks). We must note here that the 
excess of dislocations of a certain chirality leads in smec- 
tics to the break up of the homogeneous lamellar phase 
into domains separated by twist grain boundaries p2[ . 
We nevertheless pursue this case since it may applicable 
to the Aharonov-Bohm electron phase fluctuations in a 
type II superconductor and other systems where geomet- 
ric winding numbers play a role. 




FIG. 1: We complete the path r(t) with a straight segment 
connecting its beginning O and its end E. The winding angle 
fluctuation along this closed path can be calculated by de- 
composing the path into a union of non self-intersecting loops 
(denoted by the solid, dashed and dotted lines). 

When the densities of the positive and negative dislo- 
cations are equal, a more subtle averaging must be per- 
formed since the average dislocation charge in a closed 
loop is now zero. We must compute the variance of the 
dislocation charge within a closed loop. As we discuss 
below, this variance is proportional to the unsigned area 
of the loop which has a simple geometrical interpretation. 

It turns out that when the dislocations are thought 
of as point objects, it is not possible to average the un- 
signed area over random walks. The variance of the layer 
number obtained in this fashion is logarithmically diver- 
gent. Divergences are common in the statistics of winding 



numbers of random walks |l3|, |TJ, [l5| For example, 
the dispersion of the winding number of a random walker 
around a point is divergent if the walk is continuous. This 
divergence arises due to the contributions to the wind- 
ing number from trajectories which wind tightly around 
the point. The nature of our divergence is similar. By 
traversing a short distance around a dislocation a tracer 
particle travels far in the direction normal to the layers, 
and that leads to the anomalously fast diffusion in the 
direction perpendicular to the layers. 

We regularize this divergence by noting that the core 
size a determines the distance of the closest approach of 
the tracer to the dislocations. Therefore small loops in 
the trajectory are irrelevant for our purposes. We should 
therefore look at an effective discrete random walk whose 
steps are of length a taken every a 2 /D seconds. We then 
are able to calculate the variance of the layer number 
which grows as a(R, T) oc TlogT. Since a(R, T) is inde- 
pendent of i?,, the shape of a spreading cloud in this case 
is an ellipsoid which elongates in the direction normal to 
the layers. Note also that this divergence would have ap- 
peared in the case of different densities of positive and 
negative dislocations. It leads to a correction of order 
TlogT. 

We now describe our methods and results in more de- 
tail. Let r a (t) be the Brownian trajectory of the tracer 
(here a is the two dimensional vector index) . We take its 
velocity r a (t) to be random and white noise correlated 
in time, neglecting possible correlations on time scales of 
the order of the scattering time of the walker, which is 
much smaller than all other time scales in our problems. 
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Eq. (0) involves constant terms in addition to the stan- 
dard S(t — t') one. This is because averaging in Eq. (|l|) 
is done with the boundary conditions r(0) = 0, r a (T) — 
R a , needed to compute the layer number n as a function 
of position R and time T. 

Dislocations of charge qi are located at x % a . takes on 
values ±1. The layer number n(R,T) can be expressed 
in the following way 



dt 
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where e a p is the antisymmetric tensor of rank 2. Indeed, 
the expression to be integrated over time in Eq. (|^) is 
just the sum over all the dislocations of the infinitesi- 
mal change of the angle between the x-axis and the vec- 
tor connecting the tracer particle and the dislocation. 
Thus (^) is the cumulative winding number of the tracer 
around the dislocations. According to the preceding dis- 
cussion (see Fig. |l|), the function r(t) in (||) consists of 
two segments. The first is a Brownian walk from t = 
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to time T. The second is a straight line from r(T) = R 
back to the origin. 

Since we are interested in the statistical properties of 
n(R, T), expression (||) must be first averaged over dislo- 
cation strengths qi and positions xi and then over Brow- 
nian trajectories r(t). To perform the averaging, we 
assume that positive and negative dislocations are dis- 
tributed uniformly. If the total density is / = /+ + /_, 
the dislocation strength is qt = 1 with probability /+//, 
and qi = — 1 with probability /-//. 

Averaging (^|) over the strengths and positions of the 
dislocations we arrive at 

n{R, T) = (/+-/_) J dt € -f- r a f p . (3) 

The integral in (^) can be interpreted geometrically as 
the overall area covered by a vector connecting the tracer 
particle to the origin as the particle moves along its tra- 
jectory. The area is computed with the sign, so that when 
the vector rotates clockwise, the area it covers is added, 
while when it moves counterclockwise, it is subtracted 
from the answer. We refer to the integral in (||) as the 
signed area. 

Eq. (g) has a simple intuitive interpretation. If the 
densities of the positive and negative dislocations were 
the same, we would expect n(R,T) to vanish, because 
on average the tracer would encircle an equal number 
of positive and negative dislocations. The signed area 
times the difference in the dislocation densities is on av- 
erage precisely the overall number of dislocations encir- 
cled clockwise minus the overall number of dislocations 
encircled counterclockwise, which should give n(R,T). 



FIG. 2: Five equidistant in time snapshots of the isodensity 
line of the vertical slice through the expanding tracer cloud 
when /_ / /+. The cloud's shape is a figure of revolution 
of this slice around the z-axis. The lines are drawn at the 
level where the density is equal 0.1 of the maximum density 
in the center of the cloud. Note that that cloud develops into 
a biconcave shape elongated in the vertical direction (normal 
to the layers). 

At this point, we must clearly distinguish between the 
/_ = f + and /_ 7^ /+ cases. If the difference in disloca- 
tion densities is zero, then n(R, T) computed in this way 
is zero and we must take into account the fluctuations in 
the total dislocation charge enclosed by a trajectory. Let 
us first concentrate on the case /_ ^ /+. 



With these preparations, it is now straightforward to 
average powers of n(i?, T) over the Brownian walks, with 
the help of ([!]) . The average (n) of (||) can be shown to 
be zero, owing to the straight shape of the completing 
segment OE of Fig. [I]. The standard deviation a(R, T) 
can then be computed as n 2 (R,T). However we can ne- 
glect the difference ^ n 2 \ — ( n 2 ) which can be shown to 
grow slower in time than a(R,T). We obtain 



a(R,T)^{n 2 ) = 
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This is the first of the two main results of this letter. Not- 
ing that ( R 2 ) = DT for a Brownian walker, we conclude 
that the tracer particle indeed moves superdiffusively in 
the normal direction. Furthermore, (Q) gives us a way to 
calculate the approximate shape of a tracer cloud shown 
in Fig. |. 

Moreover, it is also possible to compute the entire 
probability distribution of n(R, T) which allows us to 
determine the density of the cloud. This calculation in- 
volves averaging the exponential of over the Brown- 
ian walks using Gaussian functional integral techniques. 
The answer, given in terms of infinite products, is not 
illuminating. We only note here that the probability dis- 
tribution P(n, T) of a simpler quantity n(T), which is the 
average of n(R, T) over all positions R, can be calculated 
in closed form, 



P(n, T) 
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The situation becomes more interesting when /_ = /+ . 
In this case, to compute n 2 we need to square @ first and 
then average over positions and strengths of dislocations. 
We obtain 



n 2 {R,T) = - 
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dt I dt' r a {t)r {t') G a p(r(t)-r(t')), 

(6) 



where G a p(r) — 8 a p\og(r) ^ is often referred to as 

the 2D photon propagator. This formula represents the 
unsigned area of the loop formed by r(t). It was used 
in |l7| to compute areas formed by loops in a different 
context. 

There is an intuitive way to understand why the un- 
signed area appears in this context. It can be com- 
puted geometrically as follows. First decompose a self- 
intersecting loop into a union of non self-intersecting 
subloops (see Fig. |l|). We can then show that the variance 
of the dislocation charge enclosed by the loop is equal to 
the sum of the unsigned areas of the subloops plus the 
sum over all pairs of subloops of the areas of their inter- 
sections with a plus sign if the two subloops are traversed 
in the same direction and with the minus sign if they are 
traversed in opposite directions. 
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To simplify the task of averaging (p) over Brownian 
walks, we follow the example of Ref. jl7| and rewrite the 
photon propagator in the following equivalent way 

n?(R,T) = -£ J dt J dt'h(t) h(t')x 

*(n(t)-n(0) \r 2 (t)-r 2 (t% (7) 

The advantage of this formula over (^) is in the fact that 
r*i and r2 coordinates of the 2D Brownian walker decou- 
ple and become two independent ID Brownian walks. 

It turns out that the average of (Q) over random walks 
is logarithmically divergent at t w f'. Anticipating that, 
we only need to average (0) at t — ► if. That means, we 
can neglect all the terms in ([!]) which depend on T and 
R a , while keeping only the 6(t — if) term. We obtain 



(h(*)-r 2 (t')|> 



2D\t-t' 



(8) 



and 
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Substituting this into (0) we find the leading term 



a(R,T) = (n 2 (R,T)) = ^Jdtjdt'-^ V] . (10) 

It is clear that the t = t' divergence in © should be 
regularized to give 



a(R,T) 



fDT (T 



(11) 



The regulator e appears due to the fact that the trajec- 
tory cannot wind around a given dislocation tighter than 
the dislocation core size a. The continuous formula (||) 
breaks down at distances smaller than a. The length a 
corresponds to time interval e = a 2 / D which is the av- 
erage time it takes the random walker to diffuse across a 
dislocation core. We obtain 



a(R,T) = 



fDT fDT 
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This is the second main result of this letter. Notice that 
a(R, T) does not depend on R, which should be con- 
trasted with (^) which is valid when / + =/= /_. 

To summarize, we considered tracer diffusion in a lay- 
ered system with screw dislocations. When the trans- 
verse diffusion coefficient of the is small compared to the 
in-plane diffusion coefficient, tracer particles are trans- 
ported in the direction normal to the layers by encircling 
screw dislocations. We predict the shape of a cloud of 



the tracer particles as a function of time. We find that 
size of the cloud in the direction normal to the layers (its 
height) grows faster than its width. 

To make quantitative predictions we need to address 
the following concern. The conventional transverse tracer 
diffusion coefficient D± is never identically zero. Tracer 
particles can be transported along dislocation cores or 
point defects such as pores, necks and passages as sug- 
gested by Constantin and Oswald in Q . They measured 
transverse diffusion in a thin sample of lamellar phase of a 
surfactant/ water mix. Since their sample contained only 
a few dislocations across its thickness, our effect would 
therefore not be operative. Instead one can estimate the 
effect of an isolated screw dislocation using the classic 
result of the statistics of winding numbers (e.g. to 
be negligible compared with conventional diffusion D±. 

Going back to a layered system with many screw dis- 
locations, we need to estimate the time after which the 
superdiffusion due to dislocations will dominate conven- 
tional transverse diffusion. We consider the case of equal 
densities of positive and negative dislocations because un- 
less there is a process at work which selects dislocations 
of a certain charge, the difference in the densities will be 
small. The height of the cloud due to conventional diffu- 
sion is roughly equal to the height due to superdiffusion 
when 



D±T ~ d 2 fDT In 



DT 



(13) 



where d is the interlayer distance. Assuming that the 
dislocation core size is equal to the interlayer spacing we 
obtain the crossover time 



T c 



D 



exp 



d 2 f 



(14) 



If this time is comparable to the experimentally available 
time, our phenomenon should be observable. 

Ref. M measured the anisotropy of the diffusion co- 
efficient in clean samples of diblock copolymer to be 
D±/D i=s 10 -2 . Therefore, in order for our effect to 
manifest itself, the defect density must be two orders of 
magnitude larger than d 2 f « 10~ 5 observed in the shear 
aligned diblock copolymer system of Ref. Q . The diffu- 
sion of water mixed with egg phosphatidylcholine || is 
even more anisotropic D±/D w 10~ 3 so that our effect 
can be observed for smaller defect densities. A promising 
system is a mixture of lipid and surfactant which under- 
goes a lamellar to nematic transition via proliferation of 
screw dislocations 

To conclude we mention two phenomena which require 
a modification of our predictions. First, if screw dislo- 
cations are mobile, they will tend to form bound dipole 
pairs of size comparable to the core size. Since the bound 
pairs do not contribute to the transverse transport of the 
tracer, only the density of free dislocations must be used 
in Eq. (O) . In addition, the dislocation motion 1TTI, Il9|] 
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will lead to an additional mechanism for normal trans- 
port of the tracer. Second, the presence of edge disloca- 
tions impedes in-plane diffusion of the tracer. This fact 
may be successfully taken into account by renormalizing 
the in-plane diffusivity. 

The authors are grateful to R. Sclinger for seeding idea 
which led to this calculation. 
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